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1. Introduction 

Weak decays of heavy hadrons directly probe flavor changing interactions. In the 
standard model of electroweak interactions these decays provide a means of extracting 
fundamental parameters of the theory: four mixing angles, \Vcd\, \Vcs\, \Vub\ and |Vcb|, 
which determine the decay rates of the heavy hadrons. In practice, however, our inability 
to compute these rates from first principles limits efforts to extract these fundamental 
angles. Even when limiting our attention to semileptonic decays, the simpler hadronic 
matrix elements prove too difficult to compute. 

In processes involving the semileptonic decay of a heavy quark to another heavy quark, 
as in 6 ^ ciu, a systematic expansion in powers of the hadronic scale over the heavy masses 
allows computation of the hadronic matrix elements involved at one kinematic end-point, 
namely, the point at which the resulting hadronic system is at rest in the rest frame of 
the original hadron. This "heavy quark effective theory" (HQET) is therefore of great 
relevance to the determination of | Vet \ ■ 

But for decays of a heavy quark into a light one, there is unfortunately little that the 
HQET has to say. While one can find relations between different measurable decays, it 
does not seem possible to calculate the matrix elements themselves. 

In view of the lack of first-principles calculations of heavy-to-light decay rates, phe- 
nomenological estimates depend on hadronic models used. It is no surprise that the results 
vary considerably. 

Present hadronic model calculations of heavy-to-light decay rates attempt first to 
estimate the hadronic matrix element at one kinematic point and then extrapolate based 
on some assumed functional shape of the form factor, chosen in a somewhat ad-hoc manner. 
This 'shape' is one of three types: 

(i) Pole- Dominated shapes assume the form, although not necessarily the residue, of a 
form factor dominated by the lowest- lying resonance that couples to the weak current, 
e.g., vector- meson dominance of the B* in b —>■ u transitions[|l|-|^. 
(n) Quark- Model shapes are suggested by an extrapolation of the non-relativistic con- 
stituent quark model from the low recoil region to the relativistic region. The extrap- 
olation is not unique since the non-relativistic computation does not automatically 
yield a Lorentz-invariant form factor [p|-|TT|. Relativistic wavefunctions have also been 
used to compute the shape [1T2|-|T6|. 
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(iii) Two-Component models assume a linear combination of the two preceding shapes [|T7|- 

m. 



From the practical point of view this situation is quite inadequate; from a theoretical 
one it is far worse: it leaves us with no insight into the dominant structures and dynamics 
of these hadrons. More optimistically, detailed knowledge of the correct q^ dependence 
promises both to help in the experimental determination of fundamental parameters and 
to bridge the variety of theoretical ideas underlying the present matrix element calculations. 

In this paper we argue that there is compelling evidence in favor of the first, pole- 
dominated, structure for this form factor. We show that the q^ dependence is of the 
single-pole form in the combined heavy quark, large Nc and chiral limit, where the hadronic 
matrix element is exactly computed, as we have discussed in a recent letter ETI]. We go on 



to study decays of heavy mesons in two-dimensional planar QCD, where we find striking 
confirmation of both the general results and the expected approach to the limiting case. 
While the 't Hooft model is in no sense a controlled approximation to the four-dimensional 
world, it offers valuable insights. Most importantly it provides a testing grounds for new 
models since, as we explain, many considerations motivating the modeling of form factors 
are independent of those details (dimension and number of colors) which make the model 
solvable. 

We have also discovered, in the course of this work, a number of useful relations that 
permit calculation of some matrix elements in this model for the first time. Thus we expect 
these results to be interesting as well to those interested in two-dimensional field theory 
for its own sake. 

The paper is organized as follows. In section 2 we discuss physical arguments that 
suggest the validity of the form factor shapes described above. We pay particular attention 
to arguments due to Isgur and Wise that use simple power counting to argue in favor of 
a two-component model [p!7| , p!8[| . In section 3 we prove pole dominance at large Nc- This 
is done first for the simpler two-dimensional case for which spin is absent and the heavy 
quark limit need not be taken. The four-dimensional generalization is then presented. 
In section 4 we fieetingly review the 't Hooft model in order to introduce notations and 
conventions. In section 5 we calculate weak decay form factors and discuss the results in 
section 6. The principal field-theoretic results of the paper are contained in section 3 and 
section 5, and the impatient reader may wish to jump there directly. 



2. Models of Form Factors 

Several possible shapes for the form factors of heavy-to-light decays have been pro- 
posed in the literature. Always the idea is to extrapolate (or guess) a solid calculation of a 
hadronic matrix elements at a single kinematic point to a form factor at other momentum 
transfers. 

Pole- Dominated form factors are assumed to have the functional form 



C± 
1 - q^/iJ* 



f±i<l") = . 1...2 - (2-1) 



where fi* is the mass of the lightest state which couples to the weak current V^. The 
strength of the form factor C± is obtained from by estimating the matrix element at one 
kinematic point q^ = Q^: 

C±^f±{Q^) (i-QV^*'). (2.2) 

In S ^ TT transitions this pole belongs to the B* vector meson. (In two dimensions where 
spin is absent the corresponding state is the S-meson itself which couples directly to the 
vector current). 

Pole-dominated form factors also arise in the chiral Lagrangian approach to heavy 
quark interactions [l22|-|26|| ; however this analysis is necessarily limited to small momentum 



transfers and applies not to a large range of q^ but only to the small-recoil regime in the 
vicinity of the B* pole. 

Quark-Model forms are suggested by extrapolating a non-relativistic functional shape 
to the relativistic region. This shape is computed in the low-recoil region using the con- 
stituent quark model. The extrapolation is not unique since the non-relativistic compu- 
tation does not automatically yield a Lorentz-invariant form factor. The computation 
requires an overlap integral of the non-relativistic wave-functions for the ground state 
mesons. For example, in the rest frame of the S-meson 



(/Us + A*7r)/+ + ifJ'B - /Utt)/- = V^t^BJJ'TT / dk (p^p + k)(j)B{k) 

(/+ - /_)p = y^4;i^ldk<p:&+ k) (^^ + ^ j Mk), 



(2.3) 



where pis the spatial momentum of the 7r-meson, p= {mq/^-n-)p, (px is the non-relativistic 
wave-function of the X-meson, and the tilde marks form factors which are not Lorentz- 
invariant. By rotational invariance, f± are functions of \p\'^ only. One may construct 



Lorentz-invariant form factors by writing the momentum dependence in terms of a rea- 
sonably chosen — albeit ad-hoc — replacement |p|^ -^ 9(1'^)^ such that g{q'^) -^ |p|^ as 
q^ -^ <?max = (a*b — /^tt)^- For example, since in the S-rest-frame 

q"^ = ^l + fil- 2^bE^ , (2.4) 



where E-j^ = a/|p1^ + /U^, one may take 



^(^-(^4)^) -^l (2.5) 



or, alternatively 



g{q')^[q'-{fiB + f^.f]^ (2.6) 

It is natural to suspect that these two approaches might be quite limited and that at 
best each might be applicable to a particular kinematic regime. 

Indeed, a dispersion relation may be used to write the form factor as a sum over 
contributions of resonances 

p r 2\ JB*9ttBB* . f f 2\ t n ^\ 

f±{q ) = o_„2 +f±ii ). (2.7) 

q H'B* 

where f± is the continuum contribution,* or single poles in the narrow width approximation 
or large Nc limit: 

?/ 2\ \~^ jB*9nBn ,^ o\ 

n>0 ^ '^^ 

The pole of course dominates near the kinematic endpoint when \q^ — fi'^, \ is much 
smaller than the spacing between the B* and the next resonance that couples to the current 

By the same token, the lowest term is not expected to dominate in the opposite 
situation, when \q'^ — ^\, \ is greater than the typical spacing between resonances. In this 
case many resonances will in general contribute substantially to the form factor. This is a 
signal that one is not expanding in a useful set of degrees of freedom, and a quark model 
description may be more adequate in this case. 



* The precise nature of the continuum contribution / is not important, as we need know only 
its structure in the complex q^ plane, and the strength of the associated singularity, but not the 
precise nature of the singularity. 



Two- component models attempt to capture both these behaviors particularly when 
one is interested in a large kinematic range. The character of the actual transition is 
naturally of great interest. 

While suggestive, these simple arguments prove nothing. Pole-dominance fails, for 
example, if the lowest residue is anomalously small, so that the "nearest singularity" 
dominates only for values of | g^ — /u3^ | so tiny that they are not in the physical region. It 
could happen that the higher states, often neglected, actually have large residues which 
oscillate rapidly. Form factors of this type in fact occur in QCD; they are crucial to 
understanding transition between effective field theories and quark model physics for qq 



states of a single flavor p7|]. On the other hand the residue of the lowest lying pole could be 
so much larger than all of the rest that dominance by that pole alone could be guaranteed 
for all values of q^ < fx'^. 



Isgur and Wise argued for two-component model in heavy-to- light decays |T^. This 



has also been discussed in detail by Burdman and Donoghue [p!8| . They argue that in the 
combined chiral and heavy-quark limits the form factor is pole-dominated around g^ax — 
and only there. Now it is obvious that there is a S* pole and that this pole plays a special 
role since the physical value of l^^ax ~ A*s* I i^ ^° small. 
Consider the combined limit 

/Utt -^ 0, Mb ^ CX3, A^ = /Utt/Ub fixed (2-9) 



(the analysis of Ref. [^,[0, corresponds to A = 0). Standard power counting gives 
fs* oc M^' and gnBB* oc Mb. On the other hand, /x^, — I^^b ^ {^hY , while for higher 
resonances (n > 1) one has /U^ — /U^ oc Mb- At q'^^x — if^B — I^^-kY ^^e denominator in 
( p.7| ) scales as M^, while the denominators of f± (cf. ( |2.8| )) scale as Mb- In the combined 
limit, the form factors at g^ax ^^^ pole dominated. In fact the single pole dominates over 
a region that scales like g^ax — q"^ oc Mb- This is a small fraction of the physical region, 
which scales like M^. Thus, Isgur and Wise conclude that a two-component form factor is 
appropriate: the pole term dominates around g^ax with residue given by fB*g-KBB*, and 
the quark- model term is dominant at smaller q^ . 

Though suggestive, this scenario rests on crucial assumptions. The key question is 
over what range the single pole dominates, and this depends on dynamical calculations. 
We will shortly examine these assumptions critically in the case of a fully calculable toy 
model which satisfies the same assumptions, yet displays form factors that are truly pole- 
dominated over a much larger kinematic range. It is tempting to conjecture that this holds 
in four dimensions. Before describing our exact two-dimensional results, we explain a limit 
in which the behavior holds exactly in four as well as two dimensions. 



3. Pole dominance at large A^^ 

The matrix element can be expressed in terms of two form factors f± : 



7r 



{p')K\B{p)) = ip + p'U+{q') + ip-p'U.iq') . (3.1) 



where Vj^ = q'^vQ, q = p — p' and throughout the paper p and p' alway denote the 
momentum of the B and tt mesons, respectively. The mesons B and tc have quantum 
numbers that correspond to the valence quarks Qq and qq, respectively. 

Let us derive the pole-dominated shape of S ^ tt decay in a suitable limit, and then 
discuss how the approximations might be relaxed. 

Consider QCD in the limit of large Nc, with one heavy quark Q of mass M and a 
light quark q of mass m. The form factors in ( |3.1|) are saturated by couplings of the flavor 
changing current V^ to the Qq resonances in that channel. 

To momentarily suppress the complication of spin, we pass to two dimensions where 
we can write 

{AVAB) = ^ -^r—-2 ■ (3-2) 

For odd parity states \Bn), 

{0|K|S,)=e.A?Vn. (3.3) 

We can describe these interactions conveniently in terms an effective hadron Lagrangian, 

where a field ^Pn{x) is introduced for each meson state and Ci^t couples the mesons via 
terms 

Ant = Yl 9abc{q^)ex. dxip'^d^cp'cp'' . (3.5) 

abc 

Similarly, for even parity, 

(0|K|S„) = ?./n, (3.6) 

with couplings 

Ant = 5^ ^a6c(?') 95V V- (3.7) 

abc 
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The form factors of Eq. ( |3.1|) can then be written 

even parity 

■v2 



fng^Bn{q) >r^ fng^Bn{q){lJ^B- A) 



(3i 



/-(.^)= E ^^f!^- E 



,72 _ ,,2 
odd parity ^ *" "" even parity 






: 



which is obtained with the help of the useful two-dimensional formula 

eA.?"= [-qHp + p')x + {^il-A){p-p')An^'''PpP'a ■ (3.9) 



Note that the expansions (|3.8| ) have momentum dependent numerators, proportional to 
the off-shell three point couplings, g-KBniq'^)- Using analyticity and assuming suitable 
convergence, as will be justified later, a contour integral gives 



/+(?^)= E 



fn97TBn{Pn)K 



even parity 

,2 



fnQ-KBniiJ'n) V^ fng7TBn{lJ'n)ilJ'B - t^iz) 



(3.10) 



f-iQ')= E '22 - E 



odd parity even parity 



in which the numerators are the on-shell, physical couplings. 

We make several observations. It is large Nc which allows us to treat the resonances as 
stable without continuum couplings in Eq. (|3.2|). It selects the three point couplings to one- 
particle intermediate states. In using convergence as \q'^\ -^ 00, we make an assumption 
about the large- momentum behavior of the interactions, information unavailable from a 
low-energy analysis or standard chiral Lagrangian analysis. The shift of the numerators to 
the residues at the poles, familiar in dispersion theory, has a simple physical origin: in the 
"effective" meson field theory there is freedom to make arbitrary field redefinitions without 
changing the on-shell S'-matrix. Here that freedom is used to replace the momentum 
dependent couplings (that is, the higher derivative operators) by constants at the expense 
of shifting the coefficients of higher point functions, which in turn are down by powers 
of 1/A^. We make such a shift. Note the contrast with the use of field redefinitions by 
Georgi in Ref. ||2^, in which higher point functions were suppressed instead by powers of 
the cutoff to yield a related kind of "on-shell effective theory." 

Now in the chiral limit — the light quark mass m ^ and /U^ ^ — the decay 
constants fn remain finite while the on-shell three-point elements (ttS^J-B) vanish. This is 



a direct consequence of chiral symmetry and leads immediately to a pole-dominated form 
factor. 

To show this in detail, fix the state Bn and consider the matrix element of the light- 
light current ax = ?7a759- 

{BM\ax\B{p))=p'^Fr,ip'^) 

{0\ax\ne){neBn\B) 



E 



p'"^ - iJij 



Here the sum is over the tower of qq states, \t^i)- The form factors Fn{p''^) can be written in 
terms of the off-shell couplings gtBn{p'^) and decay constants fe of the effective Lagrangian, 
Eqs. (|3.5|) and (|377|). For parity odd i?„-states 



and for parity even states 



i^n(p'^) = y:%^ (3.13) 



Again, given suitable convergence, one may apply Cauchy's theorem to these form factors 
and replace the numerators by the corresponding on-shell expressions, i.e., p'"^ -^ /i|. 

Next, axial current conservation implies fi = unless fii = 0, so the axial current 
couples only to the massless pion, tt = ttq, and only the i = term persists in Eqs. ( |3.12| ) 
and ( |3.13| ) in the chiral limit. 

Again applying axial current conservation, 9 ■ a = 0, we have 

0=p'^F^{p'') 

f-^g-KBu if n is parity even; (3-14) 

f7T97TBni'2e^''P\Pa) if n is parity odd, 

for all three states on-shell. It immediately follows that 

girBn = 0, n 7^ 0. (3.15) 

so all the coupling constants in the effective Lagrangian (|3.5|), ( |3.7| ) except to the ground 
state vanish on-shell. The case n = is singled out because, on-shell, the factor 2e^'^p\p'^ = 
0, so it need not follow that g^BB vanishes. Indeed, it does not. 



Combining these results and introducing the couphng Qt^bb = tJ^%9-n:BB{lJ-'B) in analogy 
with the definition that is natural in four dimensions, the form factors are 

. / 2n f / 2^ fBgTrBBitJ'B)t^B fBQ-KBB I o ^ a\ 

f+{q ) = -f_{q) = -^—2 = , o (3.16) 

y PS y H-B 

This readily generalizes to four dimensions, remaining still in the limit of large Nc- The 



detailed argument for the limiting case is given in Ref. ||21|]. The vector and axial- vector 
currents are of course no longer dual, and the lightest state to couple to the ttB is the B* 
vector meson, not the B itself. The essential point is that B and B* are mass-degenerate 
in the heavy quark limit so that the light-quark axial vector current which produces pions 
can rotate the states into each other. 

As above, look in the light-light channel for a fixed state B^ and take the divergence 
in a frame where p' = 0: 

0=p {B^\ax\B) = }_^p -^—-2 = ,2 _ ,,2 (3.17) 

so that chiral symmetry again implies that {ttB^IB) = which in turn implies the vanishing 
of all g-KBn except for n = 0. The way in which this comes about is naturally different 
in four dimensions, refiecting the corresponding spins, symmetries and kinematics. The 
states Bji must be analyzed according to their spin. Only for spin one is there the possibility 
of a non- vanishing coupling pi |. 



One can write 

{0\Vx\Br,) = Uex, (3.18) 

where e^ is the polarization of the Bn and then define couplings by 

{7rBr,\B)=e\p + p')x9^^^+e\p-p')J^\ (3.19) 

Nothing is learned about the g- since B^ is on-shell. It is easy to see in the B^ rest 
frame p = p' that 

g^^^t-p = Q (3.20) 

Let us restrict attention for now to the exact chiral limit in the heavy quark (infinite) mass 
limit. Then 

e-p= \p\ cos 9, 
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where 9 is the angle between the polarization and the momentum vectors, and is generally 
non- vanishing. But from the kinematics it is also true that 

I -. ^ /^n - ^1 _ / A^/2^B* = OJI/Mq) lor n = B\ 

'^' 2iJn \An + 0{l/MQ) otherwise, ^^'^'^ 

where we introduce the mass difference A^ = iin — l^'B for n^ B* states, and take the large 
mass limit in the last equality. Therefore, taking M ^ oo, the chiral limit then yields 

(n) . - r for n = S" 



" = »™ 1^1 = 1 si" A„ otherwise. P'^^) 

So the couplings to excited states vanish, or g\^ j Qj^ -^ 0. 

What can we conclude from this? We see that in four as in two dimensions, the 
effective Lagrangian in the combined limit of M ^ cxd, m —>■ 0, and Nc —>■ oo has only a 
single coupling so that the form factors for semileptonic B ^ n decays are necessarily 
pole-dominated in all kinematic regions. Moreover, in the large A''^ limit the behavior 
as the light quark mass vanishes is expected to be smooth, so if we are not strictly at the 
chiral limit, the corrections to the pole-dominated form factor are (9(/i^/A^). 

In order to make this point abundantly explicit, we compute the form factors and 
effective Lagrangian coupling to all states below in the 't Hooft model and explicitly 
compute all terms in Eqs. ( |3.2D , (^T^l), ( |3.7D . Not only will this illustrate the proof given 
here, but it allows analysis of the corrections to the limit when m ^ 0, and study of the 
dependence on the heavy mass M. 

Mesons may be classified as heavy or light in two dimensions much as in four. In four 
dimensions the QCD Lagrangian is renormalizable and masses are considered according to 
whether they are large or small compared to the scale Aqcd- In two dimensions, QCD is 
super-renormalizable but the gauge coupling, which has dimensions of mass, plays a role 
analogous to Aqcd and serves to separate heavy from light. 

It is equally clear that these arguments no longer apply when the final state meson is 
no longer the near-massless ground state. It is of great interest, therefore, to explore the 
decays B ^ ti' , B ^f p. 
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4. The 't Hooft model 

This model has been extensively studied and our work relies on technology pioneered 
by 't Hooft[^, Callan, Coote, and Gross [pUf, and Einhornpij]. In these papers the bound 



state equations were derived; and it was shown that the scattering amplitudes — and the 
form factor in particular — can be written entirely in terms of interactions among the me- 
son bound states, with no quarks in the spectrum or in the singularity structure of the 
amplitudes. 

We recall the features of the model which make it solvable, and refer the reader to 
the original papers for details. The dynamics are defined by the Lagrangian, 

L = ^trF^ ^Y,^a{l^{id^ - g^A^) - ma)^a, (4.1) 

a 

where A^ is an SU{N) gauge field, F^^, is its field strength and ijja is a Dirac fermion of 

mass nia- In the large-A^c limit, the gauge coupling is scaled with Nc. g^ = g^N is held 

fixed as N^ -^ C)0. The label a runs over two fiavors of quark, with bare masses m and M. 

We relate this Yang-Mills Lagrangian (|4.1D to the effective meson theory hadronic 



Lagrangian ( |3.4D by the most pedestrian, concrete method imaginable: by computing the 
S'-matrix elements of physical states in the model ( [4.1| ) and identifying directly with the 
coupling functions in ( |3.4| ) which reproduce the same physics. 

The main obstacle is that there are no known analytic solutions to the bound state 
equations, and we have found neither approximation techniques nor limiting cases that 
adequately serve us in the most interesting regime. We therefore turn to numerical tech- 
niques. We compute the bound state wave functions and from them evaluate the form 
factor as well the precise couplings of the low-energy effective field theory when the quarks 
are "integrated out." The techniques we use are those developed in Ref. P7[ |. 



The theory is most conveniently quantized in light-cone gauge. Because there are no 
transverse dimensions, setting A- = eliminates the gluon self-coupling. It also serves to 
project gamma matrices onto a single component in any Feynman graph that has just gluon 
vertices and (— ) component current insertions on fermion lines. The infrared divergence 
in the gluon propagator, i//c^, is regulated by taking the principal value at the pole. 

The leading term of the l/N^ expansion is the sum of planar graphs. The equations 
for the full propagator and self-energy can be solved exactly, with an extremely simple 
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result: the net effect of all gluons starting and ending on the same fermion line is just to 
renormalize the quark mass appearing in the propagator, 

rri^ —^ m^ = m'^ - g'^/rc, (4.2) 

so the full quark propagator is 

After making this shift, all remaining gluon interactions enter as ladder- type ex- 
changes. Crossings would require either gluon self-couplings, which are absent, or non- 
planar graphs, which are higher order in l/Nc- The Yang-Mills coupling constant g has 
dimensions of mass, and we choose units such that g'^/ir = 1, leaving m'^ as the single 
dimensionless number parameterizing the theory. As is well known, there is a discrete 
spectrum of free mesons. 

$(p, q) is the full meson-quark vertex, and the wavefunction (^{p, q) is defined by 

Hp, q) = — dp+ $(p, q)S{p)S{p - q). (4.4) 

ITT J 

(j) and $ only depend on p through the variable x = p-/q-, so we denote (f){x) = (j){p, q) 
and $(x) = $(p, q)q-. In terms of (/), the bound state equation may be written as 



M^-l m^-l\ [^ dy 

X l-x ) '^ Jq {y-xY 



^iMx) = + -, Mx) - i 7 — —^My)- (4.5) 



Here (f>n{x) is a light-cone momentum space wavefunction of the n^^ eigenstate, with 
mass Hn, and x = p-/q- is the fraction of light-cone momentum carried by the heavy 
quark. The (f)n{x) vanish at the boundaries, and consistency of ( [4.5| ) requires that as 

a; ^ 0, (pn{x) ~ x^, with 

ttPm cot tiPm = 1 - M^ , (4.6) 

and similarly as x ^ 1, as dictated by the boundary behavior of the Hilbert transform. 
This equation does not have solutions in terms of known functions, but may be readily 
solved numerically. 

The range of x for the bound states is always in the interval [0, 1], and ^ = outside 
of this range; but (|4.5|) determines as well the full meson-fermion-antifermion vertex. 



dy_ 

iy-z) 



^n(^) = / TT^My), (4.7) 
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for values of z ^ [0, 1]. This includes z complex, corresponding to the general case where 
one or more of the lines of the meson-quark vertex is not on-shell in its physical region. 
$^(2;) is analytic in the complex plane, with a cut on the real axis from to 1. When 
X G [0, 1], ^n{x) is defined by the principal value prescription, and 

$„ x = -^^ + + Mx), 4.8 

V X 1 — X / 



in accordance with (|4.5[ ). Since (pnix) is finite, ^nix) has zeros where the first factor on 
the right vanishes, and these are the values x± where the quarks would be on-shell. These 
zeros of the vertex function cancel quark poles in the propagators of loop amplitudes to 
ensure that no quark singularities appear in gauge-invariant Green functions. 

All loop integrations are simplified by the fact that (pn{x) is a function of x = p-/q- 
only and is independent of p+. When wave functions and propagators appear in a loop 
integrals, only the latter depend on p_|_, so the /dp_|_ is over rational functions and can be 
computed explicitly by contour integration, leaving a single integral over one real variable. 



5. Calculating the matrix elements 

In this section we calculate the form factors f± as defined in Eq. ( |3.1| ). The current 
matrix element is given by 

^J (5.1) 

X S{k)S{k-q)S{k-p) , 

where $0 is the S-meson wavefunction, \E' denotes the pion wavefunction and $^ is the 
full current-quark-quark vertex, including all resummations of gluon exchange ||31|| . This 
vertex function $^ has an intuitive meaning in the physical channel. It can be expressed 
as a sum over a complete set of resonances which couple directly to the current. This is 
evident from the formula for the Green functions. With x — k-/q-, 



^m(^' 9) = 7/x - 7- / / _^)2 ^(y' y'^ ^^^ 



Mm 
2y'{l-y')ql 



3^-+ o„.//i ~./^^2 3y.- 



(5.2) 



and 
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G{y,y;q)=2_^^ — ^ , . ■ (5.3) 



The identification of tlie current witli tlie interpolating fields for the mesons will shortly 
because obvious. 

The asymmetry between plus and minus components, evident in Eq. (^]^), is a conse- 
quence of light-cone gauge. This has made it too difficult in previous studies to compute 
the plus component of the current directly. We have overcome this obstacle and show 
below how both plus and minus components can be used with equal ease — a necessity 
for the study of flavor- changing interactions. 

For the 'good' component of the current we have 



$_(x 









= 7- 






2 M2-1 m2-l^^^^(/,^(x) 



X 1 — X 

M2 - 1 m^-1 



X 



1 — a; 



9^-^n 



1 + E 



(5.4) 



fnHlcpnjx) 



where the last equality holds for < x < 1. The decay constant fn = /q dx(f)n{x), should 
properly be multiplied by a factor y^Nc/n which we absorb into the overall matrix element 
normalization where the A^c-dependence cancels. 

A similar expression of the 'bad' plus component of the current has not been presented 
before, presumably because of the apparent complications from the associated extra factors 



in (|5.2|) . Yet things are not as bad as they look. The bound-state equation and the parity 
relation give two important results [pl|. 



Mm 



y \^~ y I Jo 



(5.5) 



and 



G{y, y;q) = S{y-x)-(q G{x, y;q ^ 



{y-xY 

Applying these to Eq. (|5.2| ) gives 

$+ = 7++7^^(-l) 



X 



1 — X 



Jnl^r 



2qi{q^ - fJ^r 



:'^Jx] 



and for < X < 1, 



(5.6) 



(5.7a) 
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1 Mm 

= 7+ +7-77^ 



2q_ x{l — x) 

2q2 \^"^ X 1-x J ^' ' q^-^l 

Examining the final vertex expressions in Eqs. ( |5.4| ) and ( |5.7^ ), each contains a weighted 
sum of pole terms plus some inhomogeneous, non-pole terms. The pole terms are of the 
same form and are readily interpreted as the contribution from the entire tower of resonant 
intermediate states. The inhomogeneous terms, however, have no such meaning; happily, 
they cancel inside matrix elements. In $+, the 7_(_ sandwiched between propagators pre- 
cisely cancels the second term. In $_, the term unity, which appears in addition to the 
sum would generally be expected to give rise to a smooth background, quark-model type 
contribution, as it arises from a coupling of the bare current to the valence quark. However, 
it drops out of the form factors since it gives a contribution in f± proportional to 



dvM-)'^i-, )--/ dt;$o(-)^(- 



q^uj — jj^uj [I — ujJq uj 1 — u uj J^ uj 1 — u 



(5.8) 



This is zero for all q^ since the two integrals cancel each other. To see this, let v = tuj m. the 
first and let f = c<;-(-w(l— a;) in the second, and rewrite each of the $(x) = J dy (f){y)/{y—x)^ 
as a Hilbert transform. Then the first integral and minus the second are each equal to 

l-i^io Jo («-(<-1)t^)" 

Therefore the matrix element of the (-I-) component of the current involves a sum over the 
same wavefunctions as for the (— ) component, with alternating signs. 

It is now easy to combine the currents in the invariant combinations q-V and s^^q^V^,, 
which couple to states of even and odd parity, or odd and even n, respectively: 



{'K\q^V,\B) = (7r(p')k+V_ +Q_V+|S(p)) = (^| - ^|)/+(?^) + (/V-(?') 
{Aq^A^\B) = {n{p')\q+V- - q-V+\B{p)) = (2£^%K) /+(?')• 



(5.10) 



Using Eqs. (g), ( pTT^ , 

q+^-±q-^+=^--^(q'-- - -^ ^1 >; ^ili^IiTIip, (5.II) 



1 


\ 


M2- 


1 


1-x ) 


n odd ^ ^'^ 
n even 


2q- 


X 
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and applying to Eqs. ( |5.1C1| ) and ( ^.11 ) we have 

r ( 2\ _ Y^ fn9-KBn\<i ) 



even n " even n 






Here 



g-KBniQ'' 



M 



q^u - ii\/uj 



1-u 







dv(t)n{v)(t)G (-)* (^ 

UJ/ \1 — UJ 



1 



u /,,, \u)/ V 1 — u; 



and 



Ar,{q^ 



fn9iTBn{<i ) 



(5.12) 



(5.13) 



(5.14) 



The g-nBn{<f') are the invariant three-point couphngs. For q^ above threshold u = p-/q- 
is real, between and 1, 



^,^) = 1 [ 1 + ^i_i4 



1 _ 2 ( f^B + f^TT \ , I f^B l^-K 



T 



r 



(5.15) 



For even parity, 



TT 



\q'V,\B) = J2 
so the second form factor is given by 



odd 






f-iQ' 



V^ _^nif)__ ^ Anjq^) ^ 2 _ 2^ 
I odd ^ I r n ^ even ^ i r n 



where 



g-rrBniq'^) = -2/U 



1-W 



dt;(/)„(z;)(/)o(-)^ 



w/ \1 — a; 

f \ , f V — UJ 



-lljvMv)^oQi^[l 



u 



(5.16) 



(5.17) 



(5.18) 



and Bn = —fng-KBn/fJ'n fo^ ^ o^d. The (^TrBn have different kinematic prefactors for even 
and odd parity as expected from the Lorentz-invariant couplings of Eqs. ( p.7|) and (|3.5|) . 
Using Cauchy's theorem, we may write 



M"') = E t4^,. 
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St 



-q^/fi 



2/, ,2 • 

n 



(5.19) 



The needed convergence of the sum is provided by the direct analysis of the wavefunc- 



tions in Ref. [0, which can be carried immediately to the present case to give that 
/ ~ 1/1^^1^+^™, and /3m > is defined in ( [4.6|) . The right hand side is far easier to 
compute, as the residues of the poles are numbers rather than functions. What is their 
meaning? 

The residues of the poles are precisely related to the three-point coupling functions one 
needs in order to construct an effective meson Lagrangian. As in Eq. ( |3.4| ) , we introduce a 
field for each meson state and proceed to write down interaction terms which will reproduce 
the 5'-matrix elements computed from 't Hooft's Feynman rules. 

The effective substitution of the on-shell values in ( |5.19| ) corresponds to making non- 
linear field redefinitions of the meson states in order to trade higher-derivative interactions 
(which arise from Taylor expanding the 3-point functions) for higher-point interactions. 
But the higher order terms — four-point, five-point, etc. — are already suppressed by extra 
factors of Nc- We can choose fields, therefore, in such a way that the only interaction terms 
are cubic couplings with no derivatives in the large Nc limit. It is interesting to compare 
these field redefinitions with those studied by Georgi in deriving an on-shell Lagrangian 
for general low-energy effective field theories |28[ . 
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6. Numerical Results 



The method used for numerical solution of the three-point couplings was introduced 
in Ref. |27], to which we refer the reader for details. The wavefunctions were expanded in 
an appropriate Fourier series and the wave function overlaps computed numerically. The 
formulas of the preceding section do not hold for the lowest residue of each tower so these 
may be fit from sum rules along the lines of Ref. [^: For any chosen q^ > {^b + fJ-ir)^ 
, every term in Eq. ( |5.19| ), can be computed directly except for ^o(/^o)- ^'i- (|5-19| ) can 
be solved for ^o(/^o) ^^^ selecting different values of q^ provides an arbitrary number of 
checks on both the accuracy and the sum rule. We find in every case the variation to be 
in the fourth significant figure at most for q^ is chosen to be of the order of the ground 
state mass times a factor of order unity. 

We have studied two values of the light-quark mass, m = 0.56 and m = 0.1 — 
corresponding to /U^ = 3.09 and /i^ = 0.72 respectively — and numerous values of the 
heavy-quark mass M. The coupling constants QTrBn depend sensitively on the light-quark 
mass but weakly on the heavy quark mass. 

Figures 1 and 2 show the approach to the chiral limit. AnifJ-^) is plotted against 
the light quark mass m for the lowest eight resonances and it is evident that An goes to 
zero for all n j^ 0. The ground state coupling, Aq, approaches a constant. Recall that 
An = —9nBnfn/lJ'n- Siuce fn and fj,^ have negligible m dependence, the shape is identical 
to that for gnBri as well. 

These figures also show the lowest residue clearly dominating the higher resonances. 
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Figure 1. The 
residue AnifJ-n) ^s. m, the hght quark mass. The value was computed numerically for 
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m = 0.1, 0.56 and a line connecting the pairs of points drawn to guide the eye. For 
n 7^ 0, the An -^ in the chiral limit, m ^ 0. 



Xl03 



JiniPn 



0.1 




0.6 



Figure 2. The 
residue A-nifi^) ^^- ^' ^^^ light quark mass, for more resonant states, n =6—18. 

Figure 3 shows how the couphngs g-j^Bn = 9-KBn/lJ'n reach constant values independent 
of the heavy quark mass M. 

2 I , , , , , , 
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S^RR ' ^n 0.5 



itBB 



-0.5 - 
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n=4 



-1' ' ' ' ' ' ' 

10» 10' 102 103 

^ Figure 

3. Approach to the heavy quark limit: —g-KBn/l^n ^^- M for m = 0.56 f/U^ = 3.09J 
and M^ = 25, 2000, 20000, 200000. Results for m = 0.1 are similar. 

While the large coupling g-jrBB* for the B* is sufficient for it to dominate, the factors 
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of the meson decay constants, /„, further ampUfy the effect on the form factor. Figure 4 
shows the rapid faU of the fn with n. 

Finally, Figures 5 and 6 show the couplings vs. n, the state number. This again 
illustrates that the lowest pole dominates and that the pole dominance is stronger as the 
pion mass decreases, in agreement with the conclusions of section 3. 
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ical n-dependence of the decay constants /„. (Shown here for m = 0.56, M = 
20000.J 
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The B* pole dominates over a large kinematic region all the way down to q^ = 0, not 
simply the small neighborhood q^ ~ Q'^^x suggested in Refs. [0, [|18 . 



7. Conclusions 

The decay B ^ nei?e is dominated by the B* vector meson pole over a large kinematic 
range. In four dimensions, the result follows from the combined chiral, heavy-quark, large- 
Nc limit. This provides a justification for the pole-dominated shape assumed at all q"^ by 
many model calculations. In contrast with previous arguments we find no evidence of a 
quark-model regime or the need for two-component models of form factors. The reason 
for the difference may be that these arguments assumed a smooth chiral behavior in which 
the pion mass could be neglected, whereas we have seen a sensitive dependence on the 
pion mass which is in any event of order of the splitting between states and shouldn't be 
expected to be negligible compared to other relevant quantities. 

The next step is to use this pole-dominance result for extracting predictions for the 
KM mixing angle Vub- This requires four-dimensional estimates of the corrections to this 
limit. Work in this direction is in progress. If the corrections turn out to be small it may 
be possible to apply these results to semileptonic decays of D mesons, for _D ^ tt or even 
D^ K. 
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In two dimensions, we have shown that for any fixed heavy quark mass the form 
factors for B -^ neue are given by a single pole at the B. The result was derived in the 
combined large Nc and chiral limit. To study the approach to the chiral limit we wrote the 
form factors as resonant sums, gave explicit formulae for the residues in terms of overlaps 
of 't Hooft wave-functions and computed these overlaps numerically. This demonstrated 
the expected behavior with decreasing pion mass and the expected scaling at large heavy 
quark mass. 

The nature of deviations from the large A''^ limit should be more easily addressed in 
two than in four dimensions. Also, there is no reason to expect that the form factors for 
B -^ n*eue, where tc* is an excited pion resonance, should be pole dominated. Work on 
these issues is in progress. 
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FIGURE CAPTIONS 

Fig. 1. The residue AnilJ^n) ^^- "^5 ^^^ light quark mass. The value was computed numerically 
for m = 0.1, 0.56 and a line connecting the pairs of points drawn to guide the eye. For 
n ^ 0, the An ^ in the chiral limit, to ^ 0. 

Fig. 2. The residue AnifJ^n) vs. to, the light quark mass, for more resonant states, n =6-18. 

Fig. 3. Approach to the heavy quark limit: —gnBn/lJ'n '^^- M for m = 0.56 (/U^ = 3.09) and 
M2 = 25, 2000, 20000, 200000. Results for m = 0.1 are similar. 

Fig. 4. Typical n-dependence of the decay constants fn- (Shown here for to = 0.56, M^ = 
20000.) 

Fig. 5. Typical n-dependence of the pole residues AnifJ^n)- (Shown here for to = 0.56, M^ = 
20000.) 

Fig. 6. Typical n-dependence of the pole residues Ani/J^n)- (Shown here for m = 0.1, M^ = 
20000.) 
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